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SUMMARY 


This  addendum  reports  the  results  of  an  extension  to  the  study  carried 
out  at  Cornell  Aeronautical  Laboratory  for  the  U.  S.  Army  (Contract 
No.  DA  30-069-AMC-645(R))  to  determine  the  time  varying  flow  in  the 
vicinity  of  a  helicopter  rotor  in  hovering  or  forward  flight  and  having  a 
fuselage  immersed  in  the  rotor  wake.  The  purpose  of  this  extension  was 
to  compare  two  models  representing  the  fuselage.  The  first  fuselage  model, 
reported  in  Parts  1  and  2,  was  based  upon  the  assumption  that  the  fuselage 
was  immersed  in  a  constant  and  uniform  flow;  the  second  model  was  based 
upon  the  assumption  that  the  fuselage  was  immersed  in  a  constant  but  non- 
uniform  flow. 

Comparisons  between  the  two  models  are  presented. 

Also  presented  is  the  digital  program  employed  to  solve  the  problem. 


FOREWORD 


The  work  reported  herein,  performed  between  January  1966  and  April  196£, 
was  accomplished  by  the  Cornell  Aeronautical  Laboratory,  Inc,  (CAL), 
Buffalo,  New  York,  for  the  Director  of  Ballistic  Research  Laboratories, 
(BRL)  Aberdeen  Proving  Ground,  Maryland.  The  research  effort  was  per¬ 
formed  under  Contract  DA  30-069-AMC-645(R)  and  was  monitored  for  BRL 
by  Mr.  Thomas  Coyle  as  Technical  Supervisor.  Dr.  Peter  ;rimi  and 
Mr.  Andrew  R.  Trenka  of  CAL  conducted  the  study.  Mr.  Harvey  Selib 
developed  the  digital  computer  program. 

This  document  is  an  addendum  to  Parts  1  and  2  of  the  final  report  under 
the  contract.  It  describes  the  modifications  made  to  the  theory  and  the 
digital  program  to  allow  the  determination  of  the  effects  of  a  fuselage 
immersed  in  a  constant,  nonuniform  flow.  Part  1  of  the  final  report 
describes  the  development  of  the  theory,  discusses  the  results  of  the 
computation,  and  provides  a  comprehensive  discussion  of  the  work  per¬ 
formed  under  the  contract.  Part  2  of  the  final  report  describes  the  formu¬ 
lation  and  application  of  the  rotor  wake -flow  computer  program  and  is  of 
use  primarily  to  those  who  plan  to  use  the  digital  computing  program. 

CAL  Report  Numbers  have  been  assigned  as  follows: 

BB-1994-S-1,  THEORETICAL  PREDICTION  OF  THE  FLOW 
IN  THE  WAKE  OF  A  HELICOPTER  ROTOR,  Part  1  -  Develop¬ 
ment  of  Theory  and  Results  of  Computations 

BB-1994-S-2,  THEORETICAL  PREDICTION  OF  THE  FLOW 
IN  THE  WAKE  OF  A  HELICOPTER  ROTOR,  Part  2  -  Form¬ 
ulation  and  Application  of  the  Rotor  Wake  Flow  Computer 
Programs 

BB-1994-S-3,  THEORETICAL  PREDICTION  OF  THE  FLOW 
IN  THE  WAKE  OF  A  HELICOPTER  ROTOR,  Addendum  - 
Effects  Due  to  a  Fuselage  in  a  Constant,  Nonuniform  Flow 
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time  and  spacial  averaged  velocities  in  the  x  and  j.  directions 

(*,</, })  components  of  velocity  induced  by  the  fuselage 

the  (*,!/.})  components  of  velocity  induced  at  a  point  m 
by  a  source  of  unit  strength 

aircraft  weight 

rectilinear  coordinates  with  origin  in  the  tip-path  plane, 
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^-component  of  the  normalized  source  strength  per  unit 
area  of  the  element 

rotor  azimuth  angle  at  k,th  position 
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1.  INTRODUCTION 


As  was  described  in  Part  1  of  this  report  (Reference  1),  a  model  for  a 
helicopter  fuselage  was  developed  consisting  of  an  array  of  quadrilateral 
source  sheets.  This  model  formed  part  of  the  representation  used  for 
determining  the  flow  in  the  wake  of  a  helicopter  rotor.  The  strengths 
of  the  individual  source  sheets  were  assigned  by  satisfying  the  fuselage 
boundary  condition  resulting  from  a  steady  uniform  free  stream.  The 
magnitude  and  direction  of  this  free  strean.  were  determined  by  computing 
an  average  in  space  and  time  of  the  flow  obtained,  in  the  absence  of  the 
fuselage  representation,  in  the  region  the  fuselage  would  normally  occupy. 

Extensive  calculations  revealed  that  the  fuselage  was  adequately  represented 
in  this  manner  for  a  wide  range  of  flight  conditions,  provided  the  region  of 
interest  in  the  flow  was  sufficiently  removed  (on  the  order  of  one  to  two* 
tenths  of  a  rotor  radius)  from  the  fuselage.  However,  considerable  spacial 
variation  of  the  time  average  of  the  flow  over  the  region  occupied  by  the 
fuselage  was  observed  to  be  present  in  most  cases,  so  that  the  flow  was 
not  well  represented  in  the  immediate  vicinity  of  the  fuselage  surface. 

This  addendum  reports  the  results  of  work  performed  subsequent  to  com¬ 
pletion  of  the  study  reported  in  References  1  and  2.  The  objective  of  this 
continuation  was  to  develop  and  evaluate  an  improved  fuselage  represen¬ 
tation  which  takes  account  of  the  nonuniformity  of  the  time  average  of  the 
flow  imposed  on  the  fuselage.  The  procedures  used  in  implementing  the 
improved  model  are  described  in  Section  2.  The  results  of  the  computations 
performed  are  discussed  in  Section  3.  A  listing  of  the  modified  digital  pro¬ 
grams  is  presented  in  Appendixes  I  and  II.  Familiarity  with  the  develop¬ 
ments  of  References  1  and  2  is  assumed  in  these  descriptions  and  discussions. 
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2.  PROCEDURES  USED  TO  ACCOUNT  FOR 

NONUNIFORM  FLOW  OVER  THE  FUSELAGE 


To  account  for  the  nonuniformity  of  the  time -averaged  flow  over  the  fuse¬ 
lage,  it  was  necessary  to  make  changes  both  in  the  main  rotor -wake -flow 
program  and  in  the  supplemental  fuselage  program.  The  changes  were 
coded  as  alternative  procedures;  so  the  original  structures  of  the  programs 
were  retained.  The  program  changes  will  be  outlined  first,  and  their 
significance  and  implementation  will  then  be  discussed. 


Consider  first  the  revisions  to  the  supplemental  fuselage  program.  The 
supplemental  fuselage  program  which  determines  the  normalized  source 
strengths  was  modified  to  solve  for  the  source  strengths  c^,  in  a  specified 
nonuniform  flow  field.  The  set  of  equations  (Equations  (20)  of  Reference  2) 


Nf 

H  Bmn 
n.-f  n 

Nf 

E  Bmn  <7J„ 

n  •  f 

was  replaced  by  the  set 


m  -  Nf 


Nf 

&  mn 

n  -  / 


*  /,  Z,  *  *  *  Nf 


(20a) 


where 


<m>  &  <n» 

*  * 


are  known  x  and  j  components  of  the  time- 
averaged  local  velocity  at  the  m*>>  element 


are  the  direction  cosines  of  the  normal  to 
the  mfh  element  with  respect  to  the  (  z,  y  ) 
axes. 
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mn 


are  coefficients  relating  velocities  induced 
by  element  n  and  its  image  on  element  m 


Nf 


is  the  number  of  fuselage  points. 


All  other  relationships  of  the  fuselage  program  are  unchanged. 
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The  main  program  which  computed  the  position  of  the  wake  was  modified 
to  account  for  the  contribution  of  a  fuselage  in  a  steady  nonuniform  flow 
as  follows.  The  equations  (Equations  (13)  of  Reference  2) 

-  Z  v j [ -  v]  (1: 


wluch  define  x  ,  y  ,  /•  components  of  the  velocity  induced  by  the  fuselage 
at  a  point  (*»</»/)»  were  replaced  by  the  equations 


V*/*>V'i)sZI  °^(V^  * 

m*t 

m*  f 

-t. 

/r»z> 


(13a) 


where 


V*m  ’  Vim  '  VIm 


are  the  time  and  spacial  averaged  velocities 
in  the  x  and  ^  directions. 

are  the  ( *  ,  y  ,  j.  )  components  of  velocity 
induced  at  a  point  by  a  source  of  unit  strength 
at  the  element  (the  tilde  (—d  indicates  the 
velocity  induced  by  the  image  of  the  element). 


That  is,  the  quantity  [  c rK  simPly  replaced  by  as  com¬ 

puted  by  solving  Equations  (20a). 


The  modified  programs  are  utilized  in  the  following  manner.  First  a  flight 
condition  is  chosen,  by  specifying  /jl,  A,  <Xr»  etc.  Then,  the  main  program 
is  run  with  the  fuselage  representation  omitted,  until  periodicity  is  estab¬ 
lished  in  the  wake  flow.  The  computations  are  then  continued  through  one 
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period  (i.  e.  through  A/4/A/#  azimuth  positions),  computing  the  flow  at  the 
points 


,  Jm  )  >  (*m  '  ‘  f  m  .  1m  )  »  m  =  1.Z,  '  "  .  Nf  . 


It  will  be  recalled  that  these  are  the  points  and  their  images  with  respect 
to  the  fuselage  plane  of  symmetry  at  which  the  fuselage  boundary  condition 
is  satisfied.  Nonuniformity  in  the  flow  is  taken  into  account  by  computing 
a  time  average  of  the  flow  at  each  of  these  points,  and  then  requiring  the 
fuselage  source  strengths  to  be  such  as  to  just  cancel  the  component  of 
this  average  flow  which  is  normal  to  the  element.  Specifically,  the  quan¬ 
tities  and  V.tm)  of  Equations  (20a)  are  computed  according  to 


»  ^  k  )  +  (*  r*i  •  ~  y m  p  frm  , 


where  Vt  and  are  components  of  fluid  velocity  in  the  %  and  directions, 
respectively,  as  defined  by  Equations  (3)  of  Reference  2,  and  ^  denotes 
the  lcih  rotor  azimuth  position.  Note  that  the  flow  at  the  image  point  has 
been  averaged  into  V™  an<*  to  account  in  part  for  any  asymmetry  in 

the  flow  with  respect  to  the  z-$  plane. 

The  values  of  v£m)  and  V*m'*  as  so  computed  are  supplied  as  inputs  to  the 
fuselage  program  and  utilized  in  Equation  (20a).  The  right-hand  side  of 
that  equation  may  be  identified  as  the  negative  of  the  component  of  the  time- 
averaged  velocity  at  (  zm  ,  ym  ,  J„)  which  is  normal  to  element  m.  The 
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source  strengths  are  computed  by  solving  Equations  (20a)  and  are  then 
made  available  to  the  main  program  where  they  are  used  in  Equations  (13a) 
when  computing  the  flow  at  a  given  point. 
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3.  RESULTS  OF  COMPUTATIONS  WHICH  ACCOUNT  FOR 
NONUNIFORM  FLOW  OVER  THE  FUSELAGE 


GENERAL  REMARKS 

It  was  reasoned  in  Reference  1  that  the  assumption  of  uniform  flow  employed 
to  determine  the  strengths  of  the  source  sheets  representing  the  fuselage  was 
most  questionable  in  the  region  close  to  the  fuselage  when  the  rotor  was  in 
the  flight  range  O.OSs^sO.  10.  In  this  flight  range,  wake  vortices  are  swept 
directly  over  the  region  which  is  occupied  by  the  fuselage.  Thus  spacial 
variations  in  the  velocity  occur  which  are  much  larger  than  those  experienced 
at  higher  or  lower  advance  ratios.  Since  the  fuselage  influence  is  substantial 
in  this  range  of  advance  ratios,  these  large  spacial  variations  should  be 
important  in  determining  the  contribution  of  the  fuselage  to  the  net  flow. 

To  determine  the  soundness  of  this  reasoning,  the  analysis  of  two  flight 
configurations  reported  in  Reference  1,  which  had  been  based  on  a  constant 
and  uniform  flow  over  the  fuselage,  were  repeated  taking  account  of  the 
nonuniform  flow  over  the  fuselage.  The  results  of  the  computations  are 
discussed  below.  To  facilitate  comparisons  between  the  two  fuselage  models, 
the  results  are  presented  for  the  same  spacial  locations  as  in  Reference  1. 

THE  FLOW  AT  A  HIGH  ADVANCE  RATIO 

The  flight  configuration  at  an  advance  ratio  of  0.  220  was  selected  as  typical 
of  the  high  advance  ratio  flows.  This  case  was  computed  using  a  nonuniform 
flow  to  determine  the  fuselage  source  strengths.  Presented  in  Figures  1 
through  4  are  plots  of  the  spacial  distributions  of  the  velocity  components 
at  various  rotor  azimuths.  Comparisons  are  shown  between  the  velocities 
computed  using  the  uniform  and  nonuniform  fuselage  flow  models. 

As  expected,  no  appreciable  difference  was  obtained  in  the  velocities 
computed  at  distances  greater  than  0.  2R  from  the  fuselage.  The  largest 
deviations  were  obtained  in  the  Vx  and  components  in  the  immediate 
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vicinity  of  the  fuselage  (see  Figure  1).  The  computed  value  of  the  com¬ 
ponent  using  the  nonuniform  flow  model  was  found  to  be  reduced  by  approxi¬ 
mately  10%,  while  the  V  component  was  increased  by  approximately  10%. 

The  spacial  variations  in  the  velocities  are  still  primarily  attributable  to 
the  fuselage  and  the  azimuthal  variations  remained  nearly  constant. 

THE  FLOW  AT  A  LOW  ADVANCE  RATIO 

Of  the  cases  investigated  in  Reference  1,  the  flow  at  /^  =  0.  0732  was  found 
to  exhibit  a  relatively  large  dependence  upon  the  fuselage  (see  Figure  46  of 
Reference  1).  It  was  also  found  that  the  spacial  variation  of  the  downwash 
in  the  volume  occupied  by  the  fuselage  (computed  with  the  fuselage  removed 
from  the  flow)  was  large  (see  Figure  55  of  Reference  1).  Hence,  this  case 
was  selected  to  be  recomputed  using  the  nonuniform  flow  model  to  deter¬ 
mine  the  source  strengths  representing  the  fuselage.  The  results  of  the 
computations  are  presented  in  Figures  5  through  10  along  with  the  corres¬ 
ponding  results  for  the  uniform  flow  model. 

It  was  found  that  all  of  the  velocity  components  were  altered.  The  magnitudes 
of  the  differences  and  whether  the  variation  was  an  increase  or  a  decrease 
depended  upon  the  distances  from  the  fuselage  and  the  rotor  plane.  The 
velocities  computed  using  the  nonuniform  flow  model  differed  from  the 
velocities  computed  using  the  uniform  flow  model  by  8%  to  14%  in  the 
region  0. 1R  to  0.  2R  from  the  fuselage. 

The  general  character  of  the  flow  was  found,  however,  to  be  the  same  for 
both  models. 
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Figure  1  COMPARISON  OF  THE  VARIATION  OF  VELOCITY  COMPONENTS  WITH 
*  FOR  y  =  0,  >  =  -0.Y5  AND  /a  =  0.220 
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Figure  2  COMPARISON  OF  THE  VARIATION  OF  VELOCITY  COMPONENTS  WITH 
y  FOR  »  =  -0.25,  >  =  -0.45, /<  =  0.220  AND  0 
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Figure  5  COMPARISON  OF  THE  VARIATION  OF  VELOCITY  COMPONENTS  WITH 
*  FOR  y  =  0,  y  =  -0. M5  AND  p  -  0.0732 
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VELOCITY  VELOCITY  '  LOCI  TV 


Figure  7  COMPARISON  OF  THE  VARIATION  OF  VELOCITY  COMPONENTS  WITH 
y  FOR  X  =  -0.25,  £  =  -0.45,  /u  =  0.0732  AND  y  >  0 
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VELOCITY 


Figure  8  COMPARISON  OF  THE  VARIATION  OF  VELOCITY  COMPONENTS  WITH 
*  FOR  y  =  0.30,?  =  -O.W  AND  /ti  -  0.0732 
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Figure  10  COMPARISON  OF  THE  VARIATION  OF  VELOCITY  COMPONENTS  WITH 
*  FOR  tj  =  0,  y  =  -0. 15  AND  0.0732 
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4.  CONCLUSIONS  AND  RECOMMENDATIONS 


Based  on  the  comparisons  made  between  the  results  obtained  using  a  fuse¬ 
lage  model  for  uniform  and  for  nonuniform  flow,  the  following  conclusions 
were  drawn. 

No  differences  between  the  two  models  were  observed  in  the  effects  due  to 
the  fuselage  at  distances  greater  than  0.  2R  from  the  fuselage.  This  was 
true  for  both  high  and  low  advance  ratios. 

At  the  low  advance  ratio  (/ *  =  0.  0732),  the  nonuniform  flow  model  was 
found  to  yield  axial  velocities  which  were  approximately  14%  lower  at 
0.  1R  from  the  fuselage  than  the  uniform  flow  model.  At  the  high  advance 
ratio  (jj.  =  0.220),  the  axial  velocity  was  approximately  10%  lower  at  the 
same  locations. 

Except  for  small  changes  in  magnitudes,  the  character  of  the  flow  was  the 
same  for  both  models. 

The  results  of  this  study  show  that  the  effects  of  the  nonuniform  flow  over 
the  fuselage  are  significant  only  in  the  immediate  vicinity  of  the  fuselage 
(within  about  0.  1R),  Otherwise,  the  simpler  uniform  flow  model  for  the 
fuselage  is  recommended  in  view  of  the  reduction  in  digital  computation 
time. 
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APPENDIX  I 


OPERATIONAL  INFORMATION  FOR  THE  MODIFIED  MAIN  PROGRAM 


This  program  is  written  in  FORTRAN  IV,  with  the  exception  of  subroutine 
CLEAR,  which  is  written  in  MAP.  This  routine  is  used  to  initialize  storages 
to  be  zero. 

INPUTS 


CARD  1 


NB: 

NRW: 

NA: 

NPNCH: 


NOPT: 


NTAPE: 

NPRINT: 


LNCT: 

NFPT: 

NXPT: 

NPINT: 


Number  of  blades,  Ng 

Number  of  Revolutions  of  wake  per  blade,  Afc 
Number  of  azimuth  stations,  A/* 

Punch  option.  If  zero,  no  cards  are 
punched  at  the  end  of  a  run.  If  not  zero, 
all  wake  point  coordinates  and  core  sizes 
at  the  final  azimuth  position  are  punched 
on  cards. 

If  zero,  the  initial  wake  configuration  is 
computed.  If  not  zero,  initial  wake 
configuration  is  read  in. 

If  not  zero,  wake  point  coordinates  and 
velocities  are  saved  on  utility  Tape  4. 

If  NPRINT  =  1,  coordinates  and  velocities 
for  each  wake  point  are  printed;  if 
NPRINT  =  2,  those  for  every  other  point 
are  printed;  if  3,  every  third;  etc. 

Number  of  lines  desired  per  page  of  output. 
Number  of  fuselage  points,  Nf 
Number  of  points  off  the  wake  for  which 
velocities  are  to  be  calculated. 

Output  is  produced  at  intervals  of  NPINT 
steps;  i.  e.  ,  if  NPINT  =  1,  the  data  for 
each  azimuth  position  is  printed. 
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CARD  2 

NFOPT: 

Fuselage  Option 

=  0;  indicates 
v*  -  /,  Vt  =  0 

*OQ  |oo 

?•  0;  indicates  uniform 

„  flow  option. 

NOUPT: 

Output  Option 

4 

1 

=  0;  no  print-out 

i  0;  print -out  and 
punch  cards. 

CARD  3 

PSIO: 

Initial  position  of  blade  1,  ,  degrees. 

REV: 

XLAM: 

XMU: 

ALPHAT: 

Number  of  revolutions  of  rotor  for  which 

calculations  are  to  be  performed,  Nm 

A 

/«■ 

aT  (degrees) 

FACTR: 

RB: 

Factor  applied  to  t  Kf 

R/b 

CARD  4,  5, 

GAMB: 

Strengths  of  blade  1;  NA  of  them. 

Al: 

Core  sizes  at  Blade 

1 ;  (NA  of  them). 

A: 

Initial  core  sizes;  (NRW)(NA)(NB)  of  them. 

Fuselage  Data:  Four  cards  for  each  point;  (4)(NFPT)  cards  in  all. 


These  are  punched  by  the  Fuselage  Program. 


Card  1 

XBAR 

X 

YBAR 

* 

ZBAR 

i 

SIGX 

SIGZ 

Card  2 

XII 

XI 2 

XI 3 

XI4 

ETA1 

ETA2 

* 

is 

?/ 

nz 
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Card  3 

ETA3 

ETA4 

Dl 

D2 

D3 

D4 

13 

7* 

eCt 

dz 

<*3 

Card  4 

XLE 

XME 

XNE 

XLZ 

XMZ 

XNZ 

Ar 

Coordinates  of  points  off  the  wake  at  which  velocities  are  to  be  computed 
NXPT  points  in  all  (up  to  three  sets  of  coordinates  per  card): 

XIPT  Y1PT  ZIPT 

*  J/  * 

Initial  Wake  Configuration  -  Read  in  only  if  NOPT  is  not  zero. 


X:  (NRW)(NB)(NA)  of  them. 

Y:  (NRW)(NB)(NA)  of  them. 

Z:  (NRW)(NB)(NA)  of  them. 

A  listing  of  the  program  is  given  on  the  pages  which  follow. 
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tlBFTC  2ZUWV  LIST.KfT 

£_ _ CAICUIAT  OF  THC  WAKt  VORTICITY  DISTRIBUTION  ^IR  A  MCI  ICOPTFR _ 

COMMON  X  (  340  *41  ,  Y  (  340 ,4 ) , Z ( 340 , 4 )  . U  ( JA'' , 4  |  ,  V I  340, 4  )  ,  W <  140,4), 

1  GAMA  (  340,4  )  ,SFG!  340,4  )(r.<t ‘491  (  I  00)  ,A1  (  in  ))  ,  M  340,4  I  ,P|  ,RAD, 

2  VMP.XMCL  ,  XMSL  ,NR,NPW,NA,NW.NOPT  ,NTAPF , NP»I NT  ,NOVCH,PS|0, 

3  XMIJ.XLAM,  ALPHAT,PINT,PS|F,XNA.DPS!  *NHl,NWI  .  XNW  *  XNR »  T  PNN, 

,4  SAT,CA?,C1,C2,PSI  •  TP  I  ,  XI  ,T|  ,t:,<J,NPS,  JPS,  IPS,  IPS  l  ,XXX,YYV, 

s _ z  /  z  *  1st  ♦_!  no,j  flg,s  i  r,  1 ,  s  i  cz ,  s  i  G3 .  g_gg_,  oe  n  ♦  m  m  . *  ni  i  ? ,  x_nu  .  m_» 

6  X  M  X  ,  XM  Y ,  X  M  7  ,  Si  G  4  ,  ST(,  S  ,  R  R  ,  S  0  iT,  < F%  S  0  I  ,  SG ,  Rr  ,  L  P  S  ,  S  F  G  l ,  S  F G  ?  , 

7  SUH,LMCT,  XIPT(400)  ,YlPT(4nO)  ,  7|PT(4.?n)  ,VX(4nn)  ,VYI 400) , 

8  VZ<400)  ,VXF(40nt  »VYF  1400)  , V/T ( 400 ) , NX PT , NAH, F ACT R , R R 
COMMON  /FUSF/  XBAR ( 100 I , YH AM ( 100 ) , Z R AR < l "0 ) , S I GX < 1 00 ) , S I GZ 1 1 00 )  , 

1  XI  1!  100)  ,XI2<  100)  ,XI  3!  100)  ,X|4(  100)  ,ETA1  (  IP**)  ,FTA?|  100)  , 

2  _ F  T  A  3  (  l  C  0  )  ,  E  TA4 1  100)  ♦  XLE  (  100)  ,  XMF  (  100)  ,XNE  (10  n|  ,XL /HOP), 

3  XXZ (  1 00  )  ,  XN Z'UOO  f, NF P T  , R  J ( 4 )  , E  J (  4 ) ,  H /( 4 )  ;'F~MJ"(  4  )  ,  D l  C  1 00  ) , 

4  0?( 100) ,D3( 100) ,041100) ,VXlNF,VYlNF,VZINF,UF,VF,WF _ 

0 IMENS I  ON  GAMRC 100) 

FQU! V ALFNCF ( GAMR 1 , GAMR ) 

DIMENSION  VXJI200) ,VZJ!?0C) 

_ 1_C  ALL  _CL  FAR1  X ,  NAB  ) _ 

CALL  CLEAR (XRAR.WF) 

CALL  CLfcARIVXJI  l),VXJ(?nOI) 

CALL  CLFARIVZJI 1  ),VZJ(200) ) 

_  PI  *  3. 1415924S36 

RAD  =  .01-74532Y25 

_ TP!  »  ?.Q*PI _ 

RF AO  1000,NR,NRW,NA,NPNCH,N0PT,NTAPE,NPRINr,LNCT,NFPT,NXPT,NP!NT, 
_  1  NDVCH  _  _  _  _ 

1000  FORMAT! 1216) 

RFAD  1000,NFOPT,NOUTP  __  _  _ _ _ 

CALL  OVOCHK ( NDVCH) 

_ IFINTAPF.L T.O)  RFW1ND  4 _ 

RFAD  100  l » PSI 0, REV, XL AM, XMU, ALPHAT , F ACTR ,RB 

1001  FORMAT!  0F8«f»)  _  ______ 

RFAD  1001, ! GAMR 1 ! I ) , I = 1 , NA ) 

READ  1001, (Al(l), 1*1, NA)  _ 

NXPT2  «  NXPT/2 

_ NB1  «  NR _ 

NW  *  NR W*NA 

NW1  *  NW*  l  _ _ _  _  _ _ 

NAR  *  NB*NA 

_ I L I M  *  NA/NR  _  _ 

l COUNT  *  0 

_ AVG  »  , 5*FL OAT! NB) /FLOAT! NA ) _ 

XNA  *  NA 

_  DPS l  *  ?.0*P1/XNA 

XNW  *  NW 
XNR  *  NR 

SAT  *  SIN! ALPHAT*RAD) 

_ CAT  «  C OS! ALPHAT*RAO) _ 

Cl  *  XMU*CAT 

C2  *  ! XMU* SAT* XL  AM) 

C3  *  XMU* SAT* SORT ! XLAM*XNB/2»0I 
XMCL  «  Cl /XL AM 
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XMSL  *  xmu*sat/xlam 

_ VX IN F  «  XM CL _ _ 

VZ  I NF  x  -FACTK*(XMSL+S0RT!,5*  XN3/XL AM)  ) 

IF  <NFOPT.NE.OI  f,0  TO  99  . .  .  _ _ _ 

VX INF  «  1.0 
V/INF  »  0.0 
99  TMP  «  Kri*DPSI 

_ TMP  I  »  SORT!  TMP*  ! TMP»  2 . 0 )  t _ _ 

FR B  *  ( T M P - T M P l Va L 0 G ( 1 . 0 ♦ T M P ♦ T M P 1 )| / 0 P S I 

READ  1 001*1  IAII  *  J )  *  I*t.NWI*J*l*NBll  .... 

IF( NFPT .EQ.O)  GO  TO  100 

REAO  1003*  !XRAP.!  I  )  *  YRAR  (  I  )  *  /BARU  I  ,  SI  GX  It)  ,  S 10/  (  I)  **LNK ,  X  IK  I) , 

1  x12m.xnm.xuu  i»etahm»fta2!II*eta3!  i i.eta4<  n* 

2  _ P I  !!_)_, J)  2  (I)  *03111  *04(1)  *XLEU)  »XMEU)  »XNE!  I)»XLZ!t  )j_ 

.3  x'xYm  ,XN2U)  ,l*l,NFPT) 

1003  F0RMATI6E12.‘>I  _  _ 

100  IF!NXPT.FQ.Ol  GO  TO  103 

READ  lOR!*(X|Prm*Y!PTm*MPTU)*t*l*NXPT)  _ 

101  on  102  I-l'NXPT 

_ CALL  FUSlGEIXlPTm*Y|PTm*ZIPTm*VXFm*VYFm*VZF<m _ 

102  CONTINUE 

_103  PS  I F  *  PS  10*  360. 0*REV _ _ _ _ 

CALL  IDOUT 

NCT  *  0  _ _ 

TPNB  -  2.0*Pl/XNB 

_ PM  »  PSI 0*RA0 _ 

PSIO  «  PSI 

PSIF  *  PSI F+R AO+O. 05  _ 

1FIN0PT.F0.0)  GO  TO  3 

2  REAO  l003»IIXII*J).l=l»NWll.J=l.NR)  _ 

READ  1003*t!Y! l*J)*I=l*NWl)  ,J*1,NBI 

_ R E_AO_JL 00 3*1  1/1  I*JI*I»1*NW1)  .Jxl.NBl _ 

GO  TO  7 

3  00  6  1=1, NWl 

XI  *  FLOAT! I-1)*0PS! 

T1  =  XI*C1 
T3  =  XI*C3 

4  00  5  J»l* NB _ _ _ 

XJ  =  FLOAT!  J-IHTPNB 

X! I  * J  )  =  COS(PSIO*XJ-XI »4T1 
Y( I  *  J I  =  SIN1PSI04XJ-XII 
Z(  I  * J  )  =  -T3 

5  CONTINUE  "  . . 

6  CONTINUE _ 

7  NPS  =  AMOO! PSI 0* 2*04 P 1 1 /OPS I  ♦  1 .5 
IF  fNPS.GT.NA)  NPS  »  1  _ 

DO  9  J* 1 *NB 

JPS  «  M00!NPS4<NA4! J- 1 1) /NB*NAB,NA) 

IF  ( JPS.FO.O)  JPS  =  NA 

IPS1  »  JPS _ 

DO  8  I«1*NW 

IPS  *  1PS1  _  _ 

ipsi  =  ips-i  "  "  .  ‘ 

IF! IPSl.EQ.O)  IPSI  «  NA  _  _ 


>• 


GAMA  (  1 1  J )  *  IGAMBI IPS1*GAM0I I  PS  l )  )  /  2 . 0 

__B_CONT  I NUF  _ 

9  CCINTINUr 

10  00  1?  J*  l  • NR  I 
00  ll  I  *  1  ,NW 

SEGl  It  J  )  =  SORT H XI  I  ,J1-X(1M  ,JI  » *  *  2*  I V  U  .  J )  -  YU  ♦  l ,  J  M  **2*I  Z I  I  .  J 1 - 
l 

11  C  ON  TJNUF _ 

12  CONTINUF 

_  13  00  29  l*l,NM _ _ _ 

DO  2fl  J*l,NHl 

XXX  =  X  (  1 1 J  )  _ _ 

YYY  *  Y ( 1 1  J  ) 

_ 7 1 1  -  Z(  I  ,  J  I _ 

lit  I*  J)  »  0.0 

V(  It  JI  =  0.0 _ _ _ 

Mil,  Jl  =  0.0 

14  00  25  L*l,NBl  _ _ _ 

1ST  *  l 

_ I  NO^NW _ 

IFLG  *  1 

IF  IL.Nf.Jl  GO  TO  16 _  _  __ _ 

INO  *  1-2 

IFLG  =2  _  _ 

IF  I  IND.G-T .0 )  GO  TO  16 

15  I S  T  =_I  ♦  l _ 

INO  =  NW 

IFLG  *  1 

IF  ( IST.GT.NW1  GO  TO  ifl 

16  SIG2  *  SORT! (XXX-XlISTtL) 1 **2* I YYY- YI IST»L1)**2*IZZ7-ZI IST»L) >  *  *  2  > 
DO  17  I R= l ST, I  NO 

_ S I G1  =  S IG2 _ 

SIG2  *  SORT  I(XXX-XIIR  +  1,L11**2M  YYV-YI  IR*l,Lll**2*l  ZZZ-Z 1 1 R*l , L 1 1 
l  **21 

SFGSO  =  SEGI I R  t L  1**2 
HMl  =  SIGl**2»SIG2**2 
IFIHMI.GT.SEGSOIGO  TO  160 

_ HN2  »  ,25*1 1SIG1»SIG2)**2-SFGSQ)*I SEGSQ- 1MG1-S1G21**?1/SEGSQ 

I F ( HM2.GT • A ( IR,L!**2 1G0  TO  160 
GGG  *  GAMA  I IR,L l/SFGl IRtLl 
GO  TO  161 

160  GGG  «  GAMAUR,L1*(  Si Gl*S IG? 1 / I S IG1 *S l G2* I  1 SIGI+SIG21 **?-SFGSQl  1 

161  XNU1  *  iYYY-YIIR*l.Lll*lZIIRtLl-ZUR*l,Ll)-IZZ7-7IIR*l,Lll* 

1 _ I  Y (  IR«L  l-YI  IR»1  »L  1 1 _ 

XNU2  *  l7ZZ-ZlIRM,Lll*CXlIR,Ll-XUR*ltlll-IXXX-XUR*ltLll* 
l  171  IR,L 1-71 IR*1 ,1  11 

XNU3  *  |XXX-XtIPM,Lll*miR,l!-YURMtLli-lYYY-YlIRM,Lll* 

l  C  X ( IRtLl-X( I R* 1 1 L 1 1 _ _ _ 

Ull,  Jl  =  Uf  I,J1*XNII1*GGG 

_ VI 1.J1  «  VI  It Jl ♦XNU2*GGG _ 

WlltJl  *  Ml l,Jl*XNU3*GGG 

17  CONTINUF 

GO  TO  118, 151, IFLG 

18  IF  (L.NE.J1  GO  TO  2> 
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IR1  =  I  -  1 

_  IF  II.FU.I)  1 R  l  -  l  _ _  —  .  .  ..  _ 

XMX  *  (Y(  IK)  ,1  |-V«  IRI^Ivl.  I  !•«/«  I«l  •  »  .»  )-;iIi*U?»i.ll-IYURl»l,ll- 

i  v  ( 1 1>  i  ♦  ? ,  i  i )  * ,  /  <  i  h  i ,  i  )  -  /  ( i « i  ♦  i .  i  n 

xmy  =  i/(ifM,u-7ii«i»i»m*ixiiKi>ifu-xnm»2,m-i;iiRi»i,u- 

1  /{|Rl*7,l  ))*<X(IKl,l>-X(tRl*l,t)> 

XMZ  *  (X(IRl,LI-X{lRlMtlH*(Y(!KUltl)-Y(lRU2,L))-(X(IRlMtU- 

_ 1  _ X(IKI»7,m*IYIIBl,l>“YIIRI»I,L)) _ 

S  IG4  *  SFG<  I R  l  ♦  1 1 L  ) 

SIG3  *  SFG< IRl.U 

SIG5  =  SQRTM  X<  l  R  l  *2 ,  U-X I  I R  1 ,  L  )  )  **?♦  IY!  I  R 1  *2 » l  )-Y(  I R 1 ,  t )  )  **2* 

1  (21 IRI»2,L)-Z(IR1.LI)**2) 

OEN  *  I  SIG3*SIG4-SIG5)*ISIG3*SI G4*S IG5) *( SI G4*$ IG5-S I G3 ) *( S IG3* 

_ I  _ SI G  5-SIG4 ) _ 

IF  (OFN.FQ.O.O)  GO  TO  25 

I F(  DEN. I  T .0 .0) WR I TE ( 6, 1002 ) l , J, St G3»  SIG4 *S IG5 
1002  FORMAT! 2X43HDEN0MINA TOR  NEGATIVE  FOR  R  COMPUTATION  I  *I3»3X3HJ  * 
l  I  3,  3X6HS IG3  *F  16. 8, 3X6HSIG4  *F  1 6. 8 , 3X6HS  IG_5  =  E16.8J^ 

RR  =  SIG3*SIG4*SIG5/SQRT( ABS(DEN) ) 

_ _ SOI l  *  SQRT ! I2.0*RR-SIG3)*1 2.0*RR»SIG3) > _ 

IF  I~S  I G  3**  2 • l E • S I G4*  *  ? ♦ S I G5**  ? )  GO  TO  19 

SF  *  !2.0»KR*S0in/SIG3 _ _ _ 

GO  TO  20 

19  SF  *  (2.0*RR-SQI11/SIG3  _  _ 

IFISF.E 0.0.0)  SF  *  L.0F-20 

20  SOI  -  SQR T 1 12.0*RR-SIG4)*! 2.0*RR»S1G4)  ) _ 

IF  I SIG4**2.LE.SIG3**?»SIG5**2)  GO  TO  21 

SG  *  ( 2.r*RR*SQI 1/SIG4 
GO  TO  27 

21  SG  =  !?.r*RR-SQl)/SIG4 
IFlSG.ro. 0.01  SG  =  l.CF-20 

22  IF  1  I  .EO.  1 )_  GO_TO_23 _ 

BF  *  (GAMA  I  I  R  1  »  J  )  *  (  A  I.  OG  {8.0*SF/A(  IR1.J)  )  ♦  .  25  )  ♦GAMAI  I  ,  J )  *  I  ALOG I  8. 0* 
l  SG/ A  (  I , J I ) +. 25)  l/(4.0*RR*S0RT!XMX**2*XMY**2«-XM7**2)  ) 

GO  TO  24 

23  RF  =  (GAMAI I ♦ J) *( ALOGf 8.0*SF/A( I,J) I ♦. 251 ) /<4.0*RR*SQRT( XMX**?» 

1  XMY**2*XMZ**2  I  ) 

24  U( I , J ]=  U( It J)»XMX»BF _ 

Vl'ltJI  =  VI  f.  J)>XM«*nF 

W( I , J)  S  W( I , J) fXM/*8F 

25  CONTINUE 

SIGl  *  SORT! XXX**2»YYY**2*ZZ7**2> 

26  00  27  1=1, NR 

_ LPS  =  MOOINPS»INA»(l-JJ  ) /N8» NA8 , NA  ) _ 

IF(LPS.FO.O)  LpS  =~  NA 
I F  f I.EO.l.ANO.L.FO.J) GO  TO  260 
PSIBK  =  FLOAT  I LPS- 1 1 *DpS I 
SINPSI  =  SIN(PSIPK) 

COSPSI  =  COSIPS  I BK ) 

RMH2  =  I XXX-COSPSI l**2*IYYY-SINpSI )**2*ZZZ**2 
IF  ( R  MM? ♦ S I G  l *  *  2  •  G  T  •  l  •  0 )  GO  10  258 
RMH  *  SORT ( RMH2  ) 

H2  =  . 25*1 l S I G1 *RMH ) **  2-1 .0 1  * ( l .0- 1  SI Gl-RMM ) **2) 

IF  (H2*RB**2.GT .1.0)  GO  TO  258 


26 


HH  *  SQRTCH2I 

_ _XHT  *  XXX»(CQSPSl»»?*SINPSI*»?/(HH»kH)  I -YYY*S I NPS I »C0SPS I ♦  (  1 « 0/ 

f  . CMHORPI-l.nl  '  '  . 

YHT  »  YYY*l SINPSI**2»C0SPSI **2/(HH*RRI 1 -X X X* S I NPS I *COSPS 1  * ( 1.0/ 
l  (HMORRI-l.ni 

7HT  a  Z///CHMOKHI 
XNIJ1  *  -YMToZ(ltLU7HT0Y(  |,L| 

XNU?  =  -7HT0XC  I  .LUXMTOZC  |  ,LI 
XND3  ®-  XHT  o  Y(  I.LIVymTOXC  1  ,L  f 

SIG2  »  SORT  ( ( XHT-X ( 1 » L  J  I ♦ *2 ♦ ( YHT- Y( 1 , l ) ) **2* C ZHT-Z C 1 , L | )*02l 
GO  TO  25 9 

250  XNU1  -  -YVYOZC 1,LI fZZ/OYC 1,U 
XNU2  »  •ZZZ*XCl(LloXXX*Z(lfLI 
XNU3  a  -XXXOYUfUoYYYOXCltH 

SIG2  »  SQKT((*XXX-X(ltLI  )**2MYYY-Y(  1*U  l**2HZZZ-ZC  I  •(.  I  >**2) 

259  GGG  *  GAMBCLPSI*<SIGUSIG2)/(SIGI*SIG2*C  <SIGl*SIG2)**2-l.0J) 
IlCItJI  =  U(  I  #  J  I  oXNUl*GGG 

\'li,  >i  -  Vv  l  ,.n 

Ml.Jl  =  WC 1 » Jl  ♦  VNU3*m,(, 

_ Gfl  TO  27 _ _ _ _ _ 

260  WI  UJI  =  W(  I»JI-GAMB(LPS)*FR8 

27  CONTINUE 

CALL  FUSLGEC  X( !  , J I , Y (I , J J , Z < I , J I , UF , VF , WF I 
UCI.JI  =  UCI , JIOXMCLOUF 
VCI.JI  *  VCI.JIOVF 

_ WCI.JJ  «  WC  I,  J|-XMSL»WF _ 

28  CONTINUE 

29  CONTINUE 

IFCNTAPE.EQ.OI  GO  TO  30 

WR I  TEC  41  PS  I » XMU, XL AH, AL PHA T , NP » NRW, N A , NW 

WRITE (A)  (CXC I. J)  ,Y(I ,JI ,Z( I, J) ,UCI ,J),V(I, JI*W(t, JI.GAMAC I, J), 
_ L  AC I.J | , I aj.NWlI f J  =  1 >  NR1 1 _ 

30  IFCNCT.NE.OIGO  TO  31 
IFCNXPT.EQ.Ol  GO  TO  3C2 
CALL  VLCTY 

IF  C  NOUTP.EU.O)  GO  TO  300 
IF  ( ICOUNT.GE. ILIM1  GO  TO  302 

300  00  3ni  I  =  1  ,NXPT2 _ 

J  =  2*1- 1 

VXJCII  a  VXJC n *AVG*(VX( JI  +  VXC jol)  I  _ _ 

VZJCII  a  VZJC IIOAVGMVZC  JIOVZC  J*l  M 

301  CONTINUE  _  _ 

ICOUNT  =  I  COUNT ♦ l 

302  C  ALL.  JOU T P IJ T _ 

31  NCT  =  NCTH 

IF (NCT.GF.NPINTINCT  =0  _  _ 

PS  I  a  PSI+DPSI 

NPS  =  NPS*  I  _ _ _ _  _ 

IF  (NPS.GT.NAI  NPS  =  1 

_ IFCPS1 .LE.PSIF I  GO  TO  3? _ 

IF(NTAPE.NE.O)  END  FILE  A 

IF(NPNCH.EQ.O)  GO  TO  320  _  _ _ 

PUNCH  1004 

1004  FORNATC  74HZZ  HELICOPTER  WAKE  VORTICITY  CALCULATIONS  -  HARVFY 
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ISEtIP  17  ) 

_  Pl»NCH_  1  ftO  I  ,  I  (  A  I  I  ,  J  )  »_!_*_1_»NW )  .  J  =  1  ,N«U _ 

PUNCH  1003.l!XU.J).I  =  l.NWll,J  =  l,NB) 

PUNCH  1^03, ((Y( I,J),I=1.NW1).J«1.NR) 

PUNCH  |  •■'0  3,  !  !Z!  I.J),  I*1.NWI )  .J=l,NB) 

320  IF  JNUUlP.tU.O)  GO  TO  1 
PUNCH  1005 

1005  FORMAT!  RCH*  ♦♦♦  *  *  »»♦»  *»♦♦»»»»»♦»»♦»♦»»♦♦♦♦♦»♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦*♦♦♦♦ 
l  ********* **'♦**'* *'ZZ_Z  1117. 1 ) 

PUNCH  1006, IVXJIII, V7JII ), 1=1, NXPT)  _  _ _ 

1006  FORMAT ( 6F 12* 5) 

WR I TF (6, 1007)  HIM, IVXJIII, VZJ1I), 1  =  1. NXPT2)  _ 

1007  FORMAT! 1HI , 39X39HVEL0C I TY  AVERAGES  FOR  OTHER  POINTS  OVFR  13. 

_ 1_ _ 9H  AZ IMUTHS//54X2HHVX,  20X2HVZ/ !  F59.5.F22.5)  ) _ 

GO  >0  l 

32  00  35  J*l,NBl _ 

TXl  =  X! 1,J) 

TYi  *  Y!  1, J) 

TZ1  *  Z(ifJ) 

_ on  33  1 =2.NW1 _ 

TX2  =  TXl 
TY?  =  TYi 
TZ2  *  TZ1 
TXl  =  X! I , J) 

TYI  =  Y! I, J) 

_ TZl  =  ZJ  T  •  J  _  _ 

X  (  I ,  J  )  =  TX2*XLAM*U( I-l.J)*OPSt 
Y  (  I.J)  =  TY?*-XLAM*V(  I-1»J)*0PSI 
Z(I,J)  =  TZ2*XLAM*W(l-l ,J)*DPSI 

33  CONTINUF 
XJ  *  FLOAT ! J- 1 ) *TPNB 

_ X (  1 , J  )  =  COS (PS1+XJ) _ 

Yll.JI  =  SIN(PSI^XJ) 

Z! 1, J)  =  0.0 

35  CONTINUF 
DO  38  J* 1 . NO  1 
JPS  =  MQO!NPS«-!NA*I  J-l)  )  7NB«-NAB,NA) 

IFtJPS.FO.O)  JPS  =  NA _ 

JPSl  =  JPS-1 

IF!  JPSl.EO.DJPSi  =  NA  _ _ _ 

SEG1  *  SFGU.J) 

GAM l  «  GAMA!  I » J I  _  _ _ 

T  A 1  *  A  (  1 ,  J  ) 

36  00  37  1=2, NW _ 

GAM2  =  GAM  1 

GAM  1  =  GAMA! I , J)  _  _ 

GAMA! I , J |  *  GAM2 
SFG2  =  SFGi 
SEG1  *  SEG! I.J) 

_ SEG(  |tJ)  =  SQRT((XU,J)-X(l»ltJ))**2»(Y(I.J)-Y!  I  ♦  1 ,  J  )  )  **2»(  Z !  I.J)- 

1  71  I ♦ 1 , J  > ) **2) 

TA2  *  TA1  _ _ _ _ 

TAl  *  A! I.J) 

A ( I.J)  =  TAT* SORT! SEG2/SEG! I.J)) 
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37  CONTINUE 

_ SEGI.1_.JI I.  jl_SORT((X_I  It  JJ-X<2.  J)  )**2MY(1,JI-Y(2.  Jl I ••?♦(*!  l.J) 

1  7(?t.lll»*2~j 

A(  1  ♦  J  I  *  AUJPS) 

GAMA(l.J)  »  (GAMB(JPS)+GAMR(JPS1)  1/2.0 

38  CONTINUE 
go  rn  13 
FNO 


U8FTC  ZZ  I OT  LIST, REF 

£ _ WAKf  VQRT1C1  TY  f.ALCULAT  1  ON  PROGRAM  -  SUBROUTINE  I  POUT _ 

SUBROUTINE  IP  OUT 

COMMON  X (140, 4 » ,Yl 340,4) ,7 1  340,41 »UI 340 ,4 ) »V 1 340,41 • W( 340,4) , 

1  GAMM34O*4|,SF(;i34O.4»,GAMBl!10P),Al  1 1 00  »  ,  A  (340 ,4  »  ,  P I ,  R  AO, 

2  VMP,  XMCL ,XMSL,NB,NRW,NA,NW,NOPT,NTAPF,NPRINT,NDVCM,PSIO, 

3  XMU,XLAM,ALPHAT,PTNT,PSIF,XNA,OPS! ,NBl , NW l , XNW, XNB , TPNB, 

4  _ S AT,r.AT,C.I  ,C?,PS1  ,TPI  ,X  I,Tl , T2 , X J  ,NPS, JPS,  IPS,  IPS  l  ,_X  X  X  ,  YYY  ^ 

5  7ZZ,  1ST,  I  NO,  IFLG,SIG1,SIG2,SIG3,GGG,DEN,XNU1  ,XNU2,XNIH,  1»1, 

6  XMX  #XMY » XMZ »  SIG4 » S IG5 , KK»  SOI l » SF  »  SQ I »  SO, BF  »L  PS » SEG1 » SFG2, 

7  SUM, INC T , XIPT(40P),YIPT( 400) »ZIPT(400)*VX(400) fVYI 400 ) , 

8  VZ( 400) ,VXF (400) , VYF ( 400) , VZF ( 400 ) , NXPT , NAB, F ACTR , RB 
COMMON  /rUSF/  XBA3  <  '  •  .  VP  AR (100 )  ,  ZBAR <100 ) ,  S I GX (1 00 )  f  s  1GZ  (100 » , 

l  XIl(100),XI.  ),Xl3(100),XI4(I00),ETAl(I0O),ETA?(I0O), 


2  E  TA3  ( 100 ) ,ETA4( 100 ) , XLF 1 100 ) , XMF ( 100 ) , XNF ( 100 ) , XLZ ( IQO ) , 

3  XX7< 100) ,XN7(100),NFPT,RJ(4) ,EJ(4),HJ(4) ,FMJ(4),01( 100), 

4  02( 100) , 03( 100),04(100),VXINF,VYINF,VZ1NF,UF,VF,WF 

1  WRITE  (A,  1000  )NR ,  NR  W,  NA ,  NFPT ,  PS  10,  PS  I F  ,  XL  AM,  XMl),  AI.PHAT,  RB,FACTR 

1000  FORMAT( IHi, 4RX33HHEL ICOPTFR  WAKE  VORTICITY  PROGRAM  //45X31HNUMBFR 

l OF  BLADES  *111  /45X31UNUMBER  OF  REVOLUTIONS  OF  WAKF  * 

2  Ill  /45X31HNUMBER  OF  AZIMUTH  STATIONS  *111/ 45X31 HNUMBpR  OF  F u S 
XELAGE  PdlNTS  *  1 1 1/45X23HPS I  (INITIAL) 

3  FU.3,8H  DEGREES  /45X23HPSI  (FINAL)  *F11.3, 

4  8H  DEGREES  /45X23HL AMBDA  =F12.5  /45X23HMU 

5  *E12.5  /45X23HALPHAT  *Fll.3,8H  DEGREES  / 

6  _  45X23HR/B _ «F 1 1 .3/45X23HVZ  INFINITY  FACTO 

7R  *F~1 1*3) 

2  DPS  *  360.0/FL0AT(NA) 

PS  *  0.0 

WRI TE( 6, 1001 ) 

1001  FORMAT ( //30X3HPS I , 15X2 2H  STRENGTH  OF  BLADE  1  ,14X20HC0RE  SIZE  AT 
l  BLADE  1  ) 


3  00  4  1*1, NA 

WR I TF ( 6, 1002  IPS, GA MB l ( I ) , A 1 ( I  )  _ 

1002  FORMAT(F35.3,E30.5,E35.5) 

PS  *  PS*DPS  _ 

4  CONTINUE 
RETURN 


END 


$ IBFTC  ZZOUTP  l  |ST,RFF 

C _  WAKF  VIIRTICITY  CALCULATION  PROGRAM  -  SUBROUTINE  OUTPUT  _ 

SUBROUTINE  OUTPUT 

COMMON  XI 34",  4  I ,Y( 140,4) , ; ( 140,4) ,01 140  *  A  I , VI 140,4 ) , W I  140,4) , 

1  GAMA  I  140,4)  ,SEG(  140 , 4  )  ,  f-AMH  1 1  ICO)  ,  A  I  (  100)  ,AI  14  0,4), PI, RAD, 

2  VMP,  XMCl.  ,XMSL,NB,NRW,NA,NW,NI)PT,NTAPF,NPR1NT  ,ndvch,psio, 

4  XMI),  XL  AM,  ALPHA?  ,  PINT,  PS  f  F,XNA,DPSI  , NB l » NWl , XNW , XN8 , TPNB, 

_ 4  _ SAT, CAT,  CJ  ,C2,PSI.TPl.Xlj>Tl,T2tXJ,NPS,JPS,  IPS,  1  PS  1 ,  XXX ,  YYY , 

5  117,  I  STVlNO,  fF'L 0,  S 1 0, f.  S l f,2  , S'fGl, "GGG, HE N,  XNU l ,  XNU2, XNU1, 1 R I , 

6  XMX,XMY,XMZ,S!G4,SIGR,RR,S0!  l,SF,SO! »SG,RF,t PS , SFGl , SFG? , 

7  SUM, LNCT , XI P T I  400 ) , Y I P T ( 400 ) , 7 1 PT I  400 ) , VX 1 400 ) , VY ( 400  I , 

8  VZ1400) ,VXF1400) ,VYF1400) , VZFI400) , NXPT , NA0 , F ACTR, RB 

1  ILINE  *  0 

_ PS  1 0  PS  I  /  RAO _ _ _ 

OP S ID  *  OPSI/RAD* 

2  DO  5  J  *  l ,  NR  I 

I  FI  II  INF.EO.O)  WR1TE(6,1000)NP,NA,NRW,XLAM,XMU,ALPHAT,DPSIO,PSID 
1000  FORMATI IHI ,46X18HHEL ICOPTER  WAKE  VORTICITY  DISTRIBUTION  //13X 

1  15HNI).  OF  BLADES  *  17,21X25HNO.  OF  AZIMUTH  STATIONS  *  II, 

2  _ 2  1X21HNU,  OF  RfV.  OF  WAKE  «  12  /9X8HL  AMBDA  Mr  12.5,  I5X4HMU  » 

3  El  2. 5, 12X9HALPHA  T  *  F7.1.5H  OFG.  UXllHOELTA  PSI  *  F7.3, 

4  5H  DEG.  //55X5HPSI  *  F8.3,8H  DEGREES  ) 

3  WRI TF I  6, 1002 1 J 

_1002  FORMATI  /59X12H8LA0E  NUMfiFR  12/1 

4  WRITE!  ft,  1.003)  (  I ,  XI  I ,  J  J  ,  Y  (  1  ,  Jl  ,  Z 1 1  ,  J 1  ,UI  I  ,  J)  ,  VI  I  ,  J )  , W1  I ,  J  ) , 

1 _ GAMA  I  1,J)  ,  A(  I  ,  J  )  ,I~l,NWl,NPRINT) _ 

1003  FORMATI  4X5HSTAT. , 10X IHX , 14X IHY,  14X IMZ , 1 4X 2HVX , 13X2HVY, 1 1X2HVZ, 

_ l  10X8HSTRENGTH,  6XRMC0RE  SIZE  /  (I  8  ,E  1 8. 5  ,JE15 . 8  I  I 

ILINE  =  11  INF  ♦  NW  l /MAXOI  NPR  l  NT  ,  l »  ♦■  3 
IF  I ILINE, GE. LNCT)  ILINE  *  0 

5  CONTINUE 

_ IFJ NXPT.FO.O)  GO  TO  6 _ 

I  FI (LINE,E0,0)WRITE(6, 1000) NB ,NA ,NRW, XL AM, XMU, ALPHAT , DPS  1 0, PS  I D 
WR!TE(A,1004)IXIPTI I > , YIPTC 1 ) , ZIPTI  l ) ,VXI I ) , VY1 I) , VZ 1 1 1 , I* l , NXPT ) 

1004  FORMATI  /83X76MVFL0CI TIES  AT  OTHER  POINTS  //19XIHX, 14XIHY, 14X1HZ, 

1  14X2HVX,13_X2HVY,J1X2HVZ  /IE26.5_,5F15._5  1J _ _ 

6  RETURN 

_ END _ 
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% IBFTC  ZZFSLG  L I  ST, REF 

C  WAKE  VORTICITY  CALCULATION  PR OGRAM  -  SUBROUTINE  FU SLGE _ 

SUBROUT  INF  FIJSLGC  ( XD » YD*  70  »Ur_D ,  VFD » WF D I 

COMMON  /FtlSF/  XBAR < 100 ) , YBAR <  100 )  , /BAR  I  1 00 ) ,  S IGXI 10O)  ,S  IGZ  I 1001  , 

1  X!  i(  1001  ,X1?I  10Q),XI3(  100)  •X|4(100)«ETA1(  100 ) «ETA2( 100) , 

2  F.TA3C  100 ) ,E  TA4(  100 1 ,  XI. F 11001  ,XME  1 1001  •  XNEI100 )  ,X|.Z  ( 1001  * 

3  XM7.(100),XN7(100),NrPT»RJ(4),EJl4),HJl4),EMJf4),OniOO)» 

_ 4 _ 07 (  100 )  « 031 100)«04(100),VXINF,VYINFtVZlNF,UF,VF,WF _ 

DIMENSION  XIK( 100*4) »ETAK( 100*41 *DDJ( 100*4) 

EQUIVALENCE  ( X I K, X II) , I ETAK ,E TA 1 ) , I DD J,D1 ) _ 

1  SUMU  *  0.0 
SUMV  «  0.0 
SUMW  «  0.0 

_ I F(NFPT.EQ.O)  GO  TO  13 _ 

2  DO  12  J* 1 » NFPT 
NFLG  «  l 

XLX  «  XME(  J)*XN7(  J)-XM7IJ)*XNEUI 
XMX  «  XNEIJ)*XL7I  JI-XN71  JI*XLFU) 

XNX  «  XLE1 J  J*XMZI  JJ-XLZI  J)*XMFU) 

_ XB* _ XO-X  BAR (J) _ 

YB  *  YD-YBARUl 

7B  *  ZO-ZRAB(  Jl  _  _ 

DS  «  (XI3(J)-X!1(J))**2META3(J)-ETA1IJ))**2 
06  *  (XI4IJJ-XI2UJI**2META4<J)-eTA2CJ>)**2 
D7  *  AMAX\(D5,06) 

3  XI  *  XLX*XB*XMX*VB»XNX*7B _ 

FTA  «  XLE( J)*XB*XMF( J ) *YB*XNE (J) *7B 

ZFTA  -  XLZI J)*XB»XMZ( J)*YB*XN7IJ)*7B 
RO  *  XI**2*ETA**2+ZETA**2 
TJ  *  R0/07 

IFITJ.LT. 6.0)  GO  TO  5 

4  S  J  (  X  I3(  J1-XII 1  J)  )»IETA2I  J)-ETA4(  J)  )  7 1  RO*SQRT(  RO ) ) _ 

VXI  *  SJ*XI 

VETA  *  SJ*ETA 
V7ETA  =  SJ^ZETA 
GO  TO  90 

5  DO  6  1=1,4 

_ R  J ( I )  =  SQRTI IXl-XlKIJtl) )**2»1ETA-ETAKIJ, I) )**2*ZETA**2) _ 

E J ( I )  =  ZETA**?MXI-XIMJ,l ))**2 

HJ1I)  «  (ETA-FTAKI J,  I))  *  I XI-XIM  J,  I  ) )_ _ 

II  =  1*1 

IFU.FQ.4)  11*1  _ _ 

TRM1  =  XI K ( J* II )-XIK|J,l ) 

_ IF J  7RM1.EQ.0.0)  TRM1  =  1.0E-6 _ 

EMJU)  =  (ETAK1 J,ll)-ETAK( J,I))7TRM1 

6  CONTINUE  _ _  _ _ _ _ 

VXI  «  0.0 

VETA  «  0.0  _ 

V7ETA  =  0.0 

7  DO  9  1  =  1,4 _ _ 

11=  I ♦ 1 

IFU.EQ.4)  11  «  1 

TRM1  «  (RJ(I)+RJ(I1  )-DDJ(  J,I))/(RJ(II«-RJ(  II)  +DDJI  J,  I )  ) 

TRM1  *  ALOG(TRMl) 
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T«M?  «  (FTAKI  J,  m-ETAM  J,m/DI)J<  J,  l> 

_ trm3  *  ( x i k i j » i )-x ik ( j  1 1 1 ) ) /nnj( J,  I  ) _ 

VX1  *  VX I ♦TRM2* TRM ] 

VFTA  *  VFTA*TRM3*TRMl 
0  IF17FTA.F0.0.0I  00  TO  9 

TRM4  *  ATAN(  <EMjm*fcJUI-HJI  m/(7FTA*RJ(  I ))) 
THMS  =  ATAN(  IEMJ(  II *F J( 11) -HJ (I l )) / ( 7ET A«R J III  II) 

V7FTA  =  VZFTA+TKM4-TRHS _ 

"  9  CONTINUF 

90  VVX  =  XlX*VX!fXlEI J)*VFTA*XL7f J)*V7FTA 
VVY  =  XMX*VXUXNF(  J)*VETA*XMZ<  J)*VZETA 
VVZ  *  XNX*VXI+XNE(  J1*VETA*XNZ(J)*V7.ETA 
GO  TO  ( 10, 11) ,NFIG 

J_0_VVVX  *  VVX _ 

VVVY  *  VVY 
VVV7  *  VVZ 
YB  =  - YD-YBAR 1 J 1 
NFLG  *  2 
GO  TO  3 

11  T R  M  =  J I G XI J)*VXINF»S 1GZ1 J)*VZ1NF _ 

SUMU'«  SUMU*TRM*IVVVX*VVX> 

SUMV  *  SUMV+TRM* ( VVVY-VVY ) _ _ _ _ _ 

SUMW  *  SU^W^TRH* ( VVVZ+VVZ) 

12  CONTINUE  _ ____ _ _ 

13  UFO  *  SUMU 

VFQ  «  SUMV _ 

WFD  *  SUMW 

RETURN  _  _ 

ENO 
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IIBFTC  Z/VLCT  t  IST.RFF 

C  WAX E  yOR T  IC  I  TV  C ALCULAT1  ON  PROGRAM  -  SUBROUT  I NF  VLCTY _ _ 

SU8ROU f INF  VI f. f Y 

COMMON  X( 340,4)  ,Y<  340,4  )  ,Z<  340,4  )  ,U|34P  ,4 1  «  VI  340.  4  I  ,  M  (140  «  4)  ♦ 
l  GAMAI340.4)  »SFG(340»4|  *GAMR1(  100)  ,  A1  (100  )  *A(  340  ,4  )  ,Pl  *RAD, 

?  VMP, XMCl ♦XMSLtNR* NR W, NA ,NW . NOPT , NT APF » NPR I  NT ,NDVCH,PSIO, 

3  XMU, XLAM, ALPHAT tPINT.PS I F.XNA, DPS! , NBl , NW1 , XNW, XNB, TPNB, 

4  _ SAT.CAT.C1.C2.PS! , TP  I , XI .T1.T2, XJ.NPS.J PS, IPS, IPS l » XXX , YYY t 

•>  ?  ZZt  1ST  *  INOt  I  FLG, SIGl , SIG2.S IG3  »GGG, DEN, XNU1  »  XNU?»  XNIJ3,  f  R  l  * 

6  XMX»XMY»XM2»S!G4»SIG5»RR»SQI1,SF»SQI» SG»  BF, L  PS, SFGl ,SEG2t__ 

7  SUM,LNCT,XlPT(4OPI,Y!PT(4O0)»ZlPTt4OO) , VXI400 ) ,VY ( 400 1 , 

8  VZ ( 400  ) , V XF ( 400 1  * VYF ! 400 1 . V Z FJ 400  K  NX P T ♦ NAB f  F ACT R * RB _ 

1  DO  7  ! *l» NXPT 

_ VXC  I  >  «  n. 0 _ _ 

VY( 1 )  *  0.0 

vzm  «  c.o  _  _ _ 

2  DO  5  J* l  »NRl 

SIG2  «  SORT!  (XIPTU  )-X(  1,J)  )**2»I  Y!  PT  t! ) -V!  1  ♦  J!)  »»?♦!  Z  IPTt  1 I- _ 

1  Z(1,J))**2) 

3  DO  4  K  =  1 j  NW _ 

SIGl  *  S1G2 

_ SIG2  *  $QRT  ( (  X I  PT( !  )-X(K»l«J)  )**2»I  Y!  PT  ID-VI  K»  1 , J  ))  »»2H  Z  IPT !  I)- 

l  7  I K+ 1  * J ) ) **21 

SFGSO  »  Sf  G( K » J I •*? 

HM1  «  S!Gl**2*SIG2**2 

_ IFJ  HM1.GT . SFGSQ ) GO  TO  30 _ 

HM2  «  .25*! ( SIGl+S I G2I ••2-SFGSOI ♦ ! SEGSO-C  SIG1~S!G2I**2)/SEGS0 
IFIHM2.GT.AIK, J)**2)G0  TO  30 
GGG  *  GAMA(K.J) /SFG(K.J) 

GO  TO  31 

30  GGG  «  GAMA(K,J)*(SlGUS!G2l/f  SIGI*S1G2*C(SIGUS!G2)**2-SEGSQII 

_ 31  X  NU  l  «  ( Y I P  TJ  IJ-Y  (K,jn*IZlK,J)-ZIK»l.  J))-(ZIPT(  M-ZIK.J)  )♦ _ 

l  (Y(K«.J)~Y(K+l«JI) 

XNU2  *  CZIPTUI-Z(K,jn*(X(Kf  J )-X t K+ 1 , J ) )•( X IPT 1 1  )-X  I K«  J )  »• 

l  (Z(Kf J I - / ( K ♦ 1 , J  >  » 

XNU3  *  (X!PTCn-X(K,JII*(Y(K,  J )-Y( K+ 1 , J )  t -(VI PTU)-Y( K, J )  t • 
l  (XIK.JI-XIKH.JI) 

_ V X (  I  )_=  VXI  n»XNUl*GGG _ 

Wi  n  *  VYI  I  )  +  XNU2*6GG 

VZ(I)  «  VZU  )*XNU3*GGG  .  _ 

4  CONTINUE 

5  CONTI NUF  _  _ 

SIGl  *  SQRTI XIPTIl)**2*YlPT(l |**2*Z!PTIl  I**2I 

_ 00  6 _ Lf  UNB _ 

IPS  *  MOIXNPSMNA*(L-in/NR*NAB,NA) 

IFILPS.EO.OI  IPS  *  NA  _ 

XNU1  *  -YIPTm*Z(l.L»*Z!PTm*Y(l,U 
XNU?  *  -MPTI!I*XII,I  dXlPTII  )*M1,U 
XNU3  «  -XIPTm*Y(l.LdYlPTm*X(l,U 

_ SIG2_L  SORT! 1X1 PT( 1 1  ~X  <  1 .1 H ♦♦?»( YI PT (  H-V(  1»L1 1»2»(2IPTC  1 )- 

l  m.LI)«*2) 

DEN  *  S!Gl*SIG2*f I S!Gl*S!G2l**2-1.0l  _  __  _ 

IF(DEN.EQ.O.O)  DEN  «  .001 

GGG  *  GAMB1CLPS»*!S1GUSIG2)/DEN 
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VX  C I )  *  VX< l )+XNUl*GGG 

_ V  Y ( I» _«  VY (!) ♦ XMU2*GGG _ 

VMM  *  VZ(  lT*XNU3*GGG 

6  CONTINUE 

vxm  «  vxi  n♦xMCL♦vxFU » 

VYIII  *  VYUUVYFt  I  ) 

VZII)  «  VZC  n-XMSL  +  VZFH) 

7  CONTINUE _ 

RETURN 

FNO 


APPENDIX  II 


OPERATIONAL  INFORMATION  FOR  THE  MODIFIED  SUPPLEMENTAL 

FUSELAGE  PROGRAM 


This  program  is  written  completely  in  FORTRAN  IV. 

INPUTS 


CARD  1 


NPTS: 

NPRNT: 

NIT: 


NOPT: 


NDVCH: 


Number  of  fuselage  elements  Nj> 

Number  of  d-j  coefficients  to  be  printed; 
i.  e.  NPRNT  =  (NPTS)  (NPTS). 

Maximum  number  of  iterations  to  be 
allowed  in  solving  the  equations  (  in  case 
of  divergence  of  the  iterations). 

^0;  indicates  nonuniform  flow  option. 

indicates  uniform  flow  option. 

Not  used. 


CARD  2,  3, 


(2)(NPTS)  +  1  yft, 

*3t  •  iil  *  ir  31  i  •  f*f  • 

xtz  •  Vtz  •  fit.  i  .  •in  •  fzt’  etc. 


A  listing  of  the  program  is  given  on  the  pages  which  follow. 
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$ IBFTC  ZZHFPF  LIST,RFF 

C  CALCULAT  l ON  Hr  POTENT (  A1  FI  CIW  ABOUT  A  HFLJCOPTCR  FlfSFI  AT.I- 

common  x i ( ino)vv~2(' YooiYoY ioo 1, <4i  fo'oT.v'i  firi'».Y?i  imi,v'3(  noi, 

1  Y4(  100 1  ,7l(  10")  ,Z2(  10"  )  ,ntP?|  (Z4(  l"0)  ,XBAR( 1"0)  ,YRAR(  100  ) 

2  «  7BAR  (  100)  ,  AMTX  C*.  4)  ,  X  PT  (  4  )  ,  YPT  (  4  )  .  7  PT  I  4  )  ,  X  1  1 1  1)0  )  .  X  !  2  (  100 ) 

3  .  X T  3 <  1001  ,X14( 100)  ,r.  i  A  1  MOO  I  ,FTA?I 100)  ,FTA31  1001  .FTA4(  10")  ♦ 

4  ZETAK  100) ,7FTA2(  100 )  ,  7F  TA3 1  1 00 )  ,7ET  A4  (  100 ) ,  XL  X  (  100  ) , 

5  _ X  MX( 100)  ,XNX(  KP)  ,*LF(  1  "_0J_,  XMF  (  ICO  )  ,XNF  (  1  "0)  ,X1  7  U  0  0  )  , _ 

6  XM7(  IOO  jVxn;  ( 10")  .R!J(  4")  ,F  fj  (4  I  ,H1  J f  4  I  ,nu  100)  ,D2(  100)  . 

7  D3( 100) .04( 100) ,F»(  100. 102)  ,SIGX(  10")  ,SIG7(  100)  ,NPTS,NDVCH, 

8  EPS, AN, BN, GN,  AX , HX , GX , AF , BF ,GE , CX ,CE ,C 7 , D5 ,0A, 07 , S J , NT  LG, 

9  FM1,EM2,FM3,EM4,XPP, YPP.ZFP, YRPP.XII J, ET AIJ , ?FTA ! J , RO, R ! 
COMMON  XfR!  J,FTARtJ,  7F:  TP  !J,  T I  J.TRIJ.VX’  ,VETA,V7ETA,  TMP1 ,  TMP?.  TMP3, 

_ L  T MP4, TMP ,VX,VY,VZ ,A1 J,AR1 J , N1 , N2 , NPRNT , MPRNT , N I T _ 

0!  MENS  I  ON  X(  100,4)  ,Y(  100,4)  ,2(  100*41  ,X  IK  (100,4)  ,FTAK(  100,4)“, 

1  7ETAK1  100,4) 

E0U1VALFNCE  ( X , X 1 ) , ( Y, Yt ) , C 7 , Z1 ) , C X IK,X 1 1 ) , C ETAK.ETA l ) , C ZFTAK, 
t  ZFTA1 ) 

OIMENSION  VXM( 100) , VZM( 100) 

1  CAL L  CLEAR!  XI, NIT) _ _ 

READ  1000, NPTS,NPRINT,N!T,NOPT, NOVCH.EPS 

l_1000  F0RMAT(5!6,F6.l)  __ 

READ  1001 , «  C  Xf !,J)*Y(!,J),Zf!,J),Jsl»4|»!*l,NPTS) 

_1001  FORMAT! 6F12.S) 

IF  (NOPT.NE.O)  READ  1001 , I VXMC I ) , VZMt!) ,1 *1 ,N«>TS) 

_ CALI  OVOCHK(NOVCH) _ 

2  DO  18  I  =  l » NPTS 

AN  *  (Y4(I)-Y2(I))*(7.3(l)-Zl(l))-(Z4(I)-Z2(l))4(Y3(I)-Yl(I)) 

RN  =  (Z4(  I)-7?(I  )  )*(X3( l )-Xl(  l )  1-1X41  I  )-X2(  I )  )♦(  731  U-/KI) ) 

GN  =  (X4(I)-X2(I>)*(Y3(l)-Yl(I))-(Y4(I)-Y2(l))4(X3( I)-Xl(I)) 
XBAR(l)  *  lXllll»X7lll»X3m*X4l!l)/4.0 

_ YBAR(l)  =  (Y1(I)»Y7(I)»Y3( l)*Y4(l))/4.0 _ 

7BAR  (  I )  *  (Zl(n  +  72(I)  +  73(I)+Z4(I))/4.0 

_ 3  IF(AN.NE.O.O)  GO  TO  8 

IE(BN.NE.O.O)  GO  TO  8 
IF(GN.NF.O.O)  GO  TO  6 

4  WR!TF(6,1002)I,(X( I , J ) , Y ( I , J ) ,2(1 *  J  )  *  J  *  1  *  4 ) 

1002  FORMAT ( 40H  BAD  SET  OF  POINTS  FOR  QUADRILATERAL  NO.  1 5/41 3F9.4 , 5X ) ) 
GO  TO  l 

6  00  7  J«l,4  _  _ 

XPT(J)  *  X(  I,J)-XBAR( I ) 

YPT( J )  »  Y( I,J»-YBAR( I » 

7PT ( J )  =  0.0 


7  CONTINUE 
GO  TO  14 


8 

AMTX( 1,1) 
AMTX( 1,2) 
AMTX ( 1,3) 
AMTX (2,1) 
AMTX (2,2) 

AN 

BN 

GN 

BN 

-AN 

• 

AMTX(2,3) 

0.0 

IF(AN.NE.O 

0)  GO  TO  11 

9 

AMTX (3,1)  « 

0.0 

AMTX( 3,21 

GN 
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AMTX (3,3)  =  -BN 

_ f)0  10  J*l,4  _  _ _ 

A«T Xf  It 4)  =  RN*YBAR(  !  )  *-GN*ZBAP (  ! » 

AMTX( ?,4J  x  HN*X( I  ,  J)  _ 

AMT  X ( 3,4)  *  GN*YU  ,  J)-BN*7( I,J) 

CALL  SIMSDL  (AMTX, 3, 3) 

XPT(J)  =  AMTXC 1 f4)-XBAR 1 1 J 

_ YPT(J)  =  AMTX( 2,4)-YBAR( T  ) _ _ 

Z  P  T I  JT  *  A  M'fx~(  3 1*4 )  -  i B  A  R  (  I) 

10  CONUNUF 
GO  TO  14 

11  AMTX  I  3.  11  *  ON 
AMTX(3,2)  =  0.0 

_ AMTX  (  3»_3 )  =_ -AN_ _ 

12  00  l  3~  jVf  ,4 

AMTX  (1,4)  =  AN*XBAR (  1  1  *BN*YRAR(  f  1  ♦GN*ZBAR  1 1 )_ 

AMTX  (  2  •  4 1  =  BN*X(ItJ1-AN*YUt  J) 

AMTX ( 3 1 4 1  *  GN*X( ! ,J)-AN*Z<  I  ,J) 

CALL  $ I MSOL ( AMT X 1 3 1  3 1 

_ X PT 1 J 1  = _ AMTX ( 1,4 1-XRAR ( I ) _ 

YPTIJI  «  AMTX( 2,4) -YBAR( I  1 

Z  PT«  J  I  =  *»MTX  (  3 ,4 1  -  7  BAR  l  1 1  _ _ _ 

13  CONTINUE 

14  AX  *  XPTI 31-XPT ( 1 )  _ 

BX  *  YPT(3)-YPT(11 

_ G X_x_Z  V T  |  3 1  - 7  PT  (JJ _ 

AF  -  HN*G X- B  X+CK 
RE  =■  C.N*AX-GX*AN 
GE  =  AN*BX-AX*BN 

CX  *  1.0/ SORT (AX**2+BX**2+GX**2)  _  _ 

CE  =  1.0/SQRT(AF**2«-BE**2*GF**2> 

_ C  2  =  1.0/SQKT( AN**2*BN*»2*GN»»2) _ 

XLXlfl  =  CXOAX 
X MX (  I  1  *  C.X*BX 

xNxn  i  =  cx*gx 
XLE1 I  )  =  CE*AF 
XME( I  )  =  CE*BF 

_ XNE ( I  )  =  CF»GE _ 

XL Z I  I  *  *  CZ*AN 
XMZI 1 1  «  CZ*BN 
XNZ( t 1  =  C7*GN 

16  DO  17  J= l , 4 

XIKIItJl  »  XLX(  I)*XPT(J|*XMX(!)*YPT(  J)  ♦XNXI  1 1  *ZPT  (  J  1 

_ E  TAK  (  I,J1  »  XLE(fl*XPT(  J)»XMFUH»YPT(  J)»XNF(  I1»ZPT(J) 

7FTAVf  I  ,  J  1  =  XLZm*XPT(  J|*XMZII1«<YPTC  JI»XNZU)*ZPT(  J1 

17  CONTINUE 

18  CONTINUE 

19  on  31  J=1 ,NPTS 

011J1  =  SORTI (XI2I JI-XI 1( J))**2*(FTA?( Jl-ETAl(Jll**21 
02(  J  )  =  S0RTnxi3(  J1-XI2I  J1  1**?*(PTA3(  J1-ETA2C  J1 1**2» 

031  J  )  *  SORT!  (XI4<J)-XI3tjnV*?META~4<J~j-ETA3(J))**?) 

04  (  J  )  x  SQRTt(XTl(J)-XI4(J))**2  +  (ETAl(J)«-ETA4(J))**?) 
05  *  IXI3I Jl-XIK J11**2+(FTA31 J1-FTA1I Jll**? 

06  r  (XI4(J)-X12(J))**?HETA4(J)-ETA2(J))**2 


38 


07  *  AMAX 1 ( D5»06 ) 

_ SJ_  *  .5»<  X  HI  Jl-Xlll  J)  )»CFTA2t JI-FTA4CJI) _ 

TRM  «  Xl?fj)-XIUJ| 

IFCTRM.E0.0.0ITRM  =  l.OF-6 _ 

EMI  *  CFTA2CJI-ETA1 C  J) l/TPM 
TRM  »  X I  3C  J  )-X  l  2 C  J  ) 

IFCTRM.FQ.O.O) TRM  =  l.OF-6 

EM2  »  CFTA3C J1-ETA2C J)) /TRM _ 

TRM  =  X I4C  J l-X 1 31 J I 
IFCTRM.FO.O.O)  TRM  *  l.OF-6 
EM3  =  CETA4C  J J-ETA3C J) ) /TRM 
TRM  *  X  1 1 ( J I  —  X 1 4 ( J ) 

IFCTRM.FO.O.O) TRM  =  l.OE-6 

_ EM 4  -  CFTA1CJ)-ETA4CJ) ) /TRM _ 

20  DO  30  I  =  l  *NPTS 

NFLG  *  1  _  _  _  __ 

XPO  *  XRARC I )-X«AR( J) 

YPP  =  YBARC  O-YBARC  J) 

7PP  =  ZBARC I I-ZBARC J) 

YRPP  *  -YBARC I )- YBAR ( J ) 

XI  U  *  XC.XC  J)*X*PP  +  XMX(  J)*YPPfXNXC  J)*ZPP 

ETA  I J  *  XLEC J)*XPP*XMFCJ)*YPP*XNF(J)*7PP _ 

7ETAIJ  *  XLZC  J)*XPP+XM7C  J)*YPP  +  XNZC'j)*7PP 

XIRIJ  *  XLXC  J)*XPP*XMXl  J)*YKPPfXNXC  JdZPP  _____ 

E1AKIJ  =  XLFC J)*XPP*XMFCJ)*YPPP*XNF I J)*7PP 

ZETRI  J  JL  XLZC  J)»YRPP+XNZC  J)»7PP _ 

Rl  *  XIRIJ**2+ETARIJ**2*-7FTRI,)**7 

RO  =  XITJ**2+ETAIJ**2+ZETAIJ**2  _  _  _ 

TIJ  =  R0/07 

TRIJ  =  RI/D7  _  _ 

21  I F { I .NE.J)  GO  TO  22 

VXI  =  0.0 _ 

VETA  =  0.0 

VZETA  *  6.2B3 1853072 
GO  TO  27 

22  IFCTIJ.GT.6.0)  GO  TO  26 

23  DO  24  K=I,4 

RIJCK)  *  SQRTCCXHJ-XIKCJ,K))»»2HFTAIJ-ETAKtJ,K))*»2»7ETAIJ»*2) 
EIJCK)  *  7ETAI J**2MXI I J-XIKC J,K) )**2 
HIJCK)  *  CETAIJ-ETAKCJ.K) )*CXl I J-XIKC J»K) ) 

24  CONTINUE 

TMPl  «  ALOGCIRlJCl)*RIJC2l-OUJ))/CRIJll)*RIJ(2)*DiCJm/OHJ) 

TMP?  *  Al0G((RtJI2)+RIJI3)-D2(jn/CRIJ(2)+RIJC3)+02(J)}>/D2(J) 

T  MP3  *  ALOGC  CRIJC3)»R1J(4)-D3CJ)  )/(-’IJC3)*R!  JC4)»D3C J)  )  )/D3C  J) 

TMP4  «  AL"oGC  CRIJC4I+RIJC  1  )-D4C J)  )  /  CR I J  C  4  I+R IJ  C  II+D4CJ) )  I/D4C  J* 

VXI  »  CETA1C  JI-ETA2C  J)  )  *TMPU  ( FT  A2 1  J ) -FT  A3  C  J  ) )  *TMP2M  ETA3C  J  )- 
l  ETA4CJ)  )*TMP3META4C  JI-ETAl  {  J)  )*TMP4 

VXI  «  -VXI 

VETA  *  CX!lIJI-X!2(jll*TMPlMX!2CJ)-X!3CJ)l  *TMP2+CX  13(J)-XI4(J)  )* 

l _ TMP3HXI4C  J)-X1I(  J)  )*TMP4 _ 

IFC7ETAI J.NE.O.O)  GO  TO  25 
VZETA  *  0.0 
GO  TO  27 

25  VZETA  *  ATANCCEMl*EIJC II-HIJC l))/C7FTAIJ*R!JCl) ) I-ATANC  CEMl*EIJ(2) 
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1  -HIJI?) |/< 7FTAIJ*RIJ!2m»ATAN<  (FM?*E!J<?1-HIJ<2I  |/<7ETAIJ 

2  *«1  J<2)»)-ATAN!  !FM2»E  I J  (  3  »-Ht  JIJH  /_!7ETAIJ*R  IJM)  ))♦ 

3  '  ATAN(  (f'H3*E  IJm-H|jl3)T/!7ETAI  j*RIJ<3MI-ATAN<  <EM3*EIJ<4) 

a  4  -HI J(4) ) /( 7FTAI J*R I J<4) ) )*ATAN< <FM4*El J<4)-HI J(4) I / I ZFTAI J 

5  *R  I  J!  4 1 )  J-ATAN!  !  EM4*E  I.M  I  1-HI  .1!  1 1 1  /  !  ZF7A I  J*R  I  J(  1 )  1 1 

GO  TO  27  _  _ _ _ 

26  TMP  =  SJ/ ! R0*SQRT  4  RO ) I 

_ VX1  *  XI I J*  TMP _ 

VFTA  *  FT  A  1‘jVt‘mp 

V7FTA  =  ZFTA I  J*TMP  _ _  _ 

27  VX  =  XLX! J|*VX1*XLF< JI*VETA*XL7(J)*V7ETA 

VV  *  XMX1 JI*VXI*XME!JI*VETA*XMZ<JI*VZETA _ _ _ 

VZ  «  XNX(J)*VXI»XN£!J1*VETA*XNZ(J)*VZETA 
_ GO  TO  <2B , 29) , NFLG _ 

28  A  I J  *  XLZ  U)*VX*XMZ!I  »*VY*XNZ  (I  »*VZ 

NFLG  *  2  _ _ _ 

X  1 1 J  *  XIRIJ 

ETAIJ  *  FTARIJ _ _  _ _ _ 

ZETAIJ  *  ZETRIJ 

_ RO_  «  R  I _ 

TIJ*  *  TRIJ 

GO  TO  2 2  _  _ _ 

29  ARIJ  *  XLZ! I  1*VX-XMZ!I )*VY*XNZI  I  )*VZ 

B ( I «  J )  *  AIJ*ARIJ _ 

30  CONTINUF  . 

31  CONTINUE _ 

N 1  *  NPTSM 

_  N2  *  NP TS ♦ 2  _ 

32  00  33  I* 1»NPTS 

IF  (NOPT.NE.Ol  GO  TO  320  _  _ 

R(  I «  N 1 )  *  -XL Z I  I  ) 

_ 9 ( I , N 2)  «  -XNZ ( I ) _ 

GO  TO  33 

320  Q  (  1  •  N I  I  *  -XlZ<l)*VXM!i  l-XNZm_*VZM!II _ __ _ 

0 ( I »  N2 )  *  B( I »N  l ) 

33  CONTINUF  _  _  4  _  _ 

IF(NPRNT.FQ.O)  GO  TO  38 

_ MPR  I  NT_  *  MINQ!NPRNT,NPTSI _ 

34  DO  37  I=l,MPR!NTfB 

IR  =  MINO! I*7,MPRINT1  _ 

WRITE<6, 1003)1, 18, MPRINT 

_1003  FORMAT!  1HI  ,44X42HP0TENTT  AL  FLOW  ABOUT  A  HJEIJ  COPTER  JO)S  El  AGE  7/5  OX 
1  11HBIJ  FOR  J  *13, 2H  -13, 9H  ,  I  «  1-13/) 

35  00  3 6 _ J  * 1 , MPR I  NT _ 

WR I TF ( 6,10041 (B(J,K)«K=l,l8)  ' 

1004  FORMAT! RE16. 5)  _ 

36  CONTINUE 

37  CONTINUE  _  _  _ _  _ _ _ _ 

38  CALL  SIMEQ 

_ C  A  LL_ OUTPUT _ 

GO  TO  1 

ENO  _ _ 
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$ IBFTC  ZZSMEQ  LIST, REF 

C  CALCULATION  OF  POTENTIAL  FLOW  ABOUT  A  HELICOPTER  FUSELAGE  -  SIHFQ 
SUBROUTINE  SIMEO 

_  COMMON_Xl(  100) ,X2( 100) , X3 ( 1 00 1 ♦ X . : ICO ) , VI C 1001 , Y2 C IPO  I f Y3( 100  I , 

l  Y4(100),ZI( 100),72( 100>,Z3(100) ,74( 130) , XDAR ( 100 ) , YBAR(IOD) 

_ 2  _  ,ZBAR(  100)  ,AMTX(3,4)  ,XPT  (4 )  ,  YPT  ( 4)  ,  ZPT(  4 )  ,  X 1  U  100 )  ,  X  1 2U00) 

3  ’  ,XI3(100),XI4(100),ETAKIOO),FTA?(100),ETA3(100),FTA4(100), 

4  _ Z_F TA1(  100)  ,7.FTA?U00)  ,  7ETA3  (  100  )  ,  ZF  TA4  (  100 ) ,  XL  X  (  100), _ 

5  XMX(iOO)  ,XNX(  100),XLF(  100)  , XMF ( 100 ) ♦ XNE (  100 )  ,  XI.Z  I  100 )  , 

_  6  XMZfiOO) ,XNZC10C).R1J(4),EIJ(4),H1J(4),0H 100 ) , 0?( 100) ,  • 

7  03(100),  04(1 00  ),B(  100,  1 02 )  ,  S 1GX  ( 100  ) ,  S IGZ  ( l  O'' )  ,NPTS,  NDVCH, 

8  EPS,AN,BN,GN,AX,BX,GX,AE.BE,GF,CX,CE,CZ,05,D6,07,SJ,NFIG*_ 

9  EMI  *  EM2*  FM3»EM4,  XPP,  YPP,  ZPP,  YRPP  »  XI  1 .1,  ETA  I  J,  ZFT  A I J  ,  RO,  R 1 
COMMON  XIRIJ,ETARIJ, ZETR I J,TI J,TR! J,VXI , VFT A, VZET A , TMPl , TMP2, THP3, 

l  TMP4, TMP , VX, VY* V7 , AI J , AR I J,Ni , N?,NPRNT, MPRNT ,N I T 

DIMENSION  X(  100,4)  ,Y(  100,4)  ,Z(  100_»VUXIK(  100,4)  ,FTAK(  100,4),  _  _ 

l  ZETAK ( 100, 4) 

EQUIVALENCE  ( X , XI ) , ( Y, Y 1 ) , ( Z , Z 1 ) , ( X  IK ,X 1 1 ) , ( ETAK, ETA  1 ) , ( ZET AK , 
l  ZETA1) 

_ 0 IMENS ION  Ul( 100) ,U2( 100) _ 

1  JS  *  Ni  v 

_ EPS l  *  100, 0*EPS  _ _ _ _ 

2  IT  =  0 

3  00  A  1*1, NPTS  _ _  _  _  _ 

Ul( I )  *  0.0 

4  00  5 _ J* 1 , NPTS _ 

TERM  =  —  B ( I , J ) / B ( 1,1) 

IFd.EQ.J1  TERM  =  0.0  _ _  _  _  _ 

Old)  *  UK  I  )*TERM*Ui(J) 

5  CONTINUE  _ _  _ _ 

U 1 ( I )  =  Ul( !)♦«( I,JS)/B(I,I ) 

6  CONTINUE _ 

7  DO  10  1=1, NPTS 

IF(U2( I ) .NE. 0.0)00  TO  8  _ 

TMP  =  AHS(U2( I ) — Ul (1)1 

GO  TO  9  _ _ _ 

8  Tmk  =  h0S( (U2( I I-U1( I) )/U2( I ) ) 

_ 9  IF(TMP.GT .EPSi  GO  to  K> _ 

10  CONTINUE 

I F f JS.EQ.N2)  GO  TO  13  _ 

11  00  12  1*1, NPTS 

S IGX (  I )  *  Ul( I )  _  _ 

12  CONTINUE 

_ JS  =  N2 _ 

GO  TO  2 

_  13  00  14  1*1, NPTS  _  _ _ 

S IGZ (  I)  *  UllII 

14  CONTINUE  _ 

RETURN 

15  IT  =  I T»  1 _ 

IF  ( (T.GE.NIT)  GO  TO  18 

16  00  17  1*1, NPTS 

U2( I )  *  Ul( I) 

17  CONTINUE 
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go  rn  3 

_ IB  IFJJS.FO.N2)  GO  TO  I **  _  _ _ _ 

WR l TF (6. 1000)  EPSl » NI T 

1000  FORMAT  ( 5X48HEQUAT IONS  FOR  SIGMA  X  010  NOT  CONVERGE  TO  WITHIN  F7.4 

1  ♦ 1 ?H  PER  CENT  IN  IA.ILH  ITERATIONS  ) 

GO  TO  11 

19  WR I TE I  6. 1001 )  EPSl. NIT 

1001  FORMAT!  SX48HEQUAT IONS  FOR  SIGMA  l  010  NOT  CONVERGE  TO  WITHIN  F7.4 

l  « 12H  PERCENT  IN  16.11H  ITERATIONS  » 

GO  TO  13  _ _ _  _ 

ENO 
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SIBFTC  7Z0UT  1. 1  ST.RFf 

c _  CALCULATION  Of  POTENTIAL  FLOW  AROUT  A  HFL ICQPTER  FUSHAGF  -  OUTPUT 

SUBROUTINE  OUTPUT 

COMMON  XI  (  IOO)  VX?(  lOO),X3(  100  1  .X4(  IOO)  ,Vl  (IDO)  (  V?M<"I,V3(  100), 

1  Y4I  100)  »7lf  100)  ,??(  10°  '  MI  100)  ,  74(  nOI  .XBARI  ICO ) , YRAR I l 00 ) 

2  ,7HAR!  IOO.AMTXI  3,4)  ,XPT 14 ) , VPT ( 4  I , ZPTI 4  I , X  I  1  I  ICO  I , XT ?I 100 ) 

3  .  XI  31100)  ,XI4I  100)  .FTAKICO)  ,rTA?I  1001  .ETA3I  |0*|  ,FTA4!  100) . 

4  _ 7 F  T A l  (  1 0 Oh  TFT A?  ( 1  0*  )  , 7.E T A 3 (  100  )  ,  7 F  T A4 (  1 00 )  ,  XLX  (  100  ) , 

5  XMXUOOV,XNX(  ICO)  ,XLF(  100 ) , XME l ICC ) , XNE I  100 » , XL 7  I 100)7 

6  XM7( 100) ,XN7I 100) ,RI J(4) ,E I J ( 4 ) ,HI JI4) ,01(10*1 ,021 100) , 

7  03(  100) ,04( 100) ,B( 100,  )C2) ,SIGXI 100) , S IG7 ( 1 00 ) ,NPTS, NOVCM, 

8  EPS, AN,BN,GN,4X,BX,GX,AE,BE.GF,CX,CE,CZ,05,06,D7,SJ,NFIG, 

9  EM1,FM2,EM3,EM4,XPP,YPP,7PP, YRPP.XI I J.ETAI J,ZETAI J,RO,RI 
_ COMMON  XI R1 J.ETARl J,Zf TR I J,T1  J,TR1  J.VXI  ,VFTA,VZETA,TMP1, TMP2, TMP3, 

l  TMP4, TMP , VX, VY, V7 ,A I J, AR1 J.Nl ,N2,NPRNT, MPRNT, N I T 

DIMENSION  X(  100,4)  ,Y(  100,4)  ,  7  (  100 ,4)  , X  I K (  1 00,4)  , F.TAK  1100 , 4 ) , 
l  ZETAKI 100,4) 

EQUIVALENCE  I  X ,  X 1 ) ,  I Y,  Y 1 )  ,  (  Z  ,  Z  1)_ ,  I X  I  K,  X  1 1 )  ,  I  EJT AK ,  FT  A  l )  , I ZETAK, 

1  ZETA1) 


1  PUNCH  1000, NPTS 

1000  EflRMAT (  32HZ 7  HELICOPTER  FUSELAGE  PROGRAM  I6.22H  POINTS  -  HARVEY 


1SEL IB, 12X2HZZ  ) 

11  *  l 

12  *  2 

_ 13  *  JL _ 

14  =  4 

2  on  3  1=1, NPTS 

PUNCH  1*01,  XBARI  I  )  ,  YBAR  (  l  ) ,ZRA«(  I )  ,  S IGX  U) ,  S IGZ  (  I  )  .ZERO,  1,11 
PUNCH  lOOl.XlK  I  ) ,  X  I2(  I  )  .  X  I  3  (  (),X(4(I  )  ,  ETA  1  (  I  ),£TA2(  I  I,  I  ,12 
PUNCH  ICCl ,  ETA3 ( I) ,FTA4 (I ) ,01 (I) , 02 (I) , D3 ( I ) , 04 ( I ) , I , I  3 
PUNCH  1 00  1 ,  XL  FI  I ) ,  XME I  I)  ,XNE(  I) ,  XLZ  II)  .  XM7  (I) ,  XNZ  II)  ,  1 ,  14 _ 

1001  FORM  ATI 6E 1 275 ,16,121 

3  CONTINUE 
NPAGF  =  NPTS/50 

IF(NPAGE*50.LT.NPTS)NPAGE  =  NPAGEH 

4  00  5  1=1, NPAGE 

_  1 1  -  JO*_l  I-  1  )  »  1 _ 

I?  =  MING (  I  1  +  43, NPTS ) 

WRITE (6, 100?) ( S IGX ( J ) • S IGZ ( J ) ,XLZ ( J) ,XMZ(J) , XNZ I J ) , J= 1 1 , 1  2) 

1002  FORMAT!  1M1,44X42HP0TENUAL  FLOW  ABOUT  A  HFL  ICOPTER  FUSELAGE  //IBX 

1  7HSIGMA  X, I 5X7HS IGMA  Z , 12X1 1HLAMB0A  7ET A ,  14X7HMU  ZFTA.15X 

2  7HNU  ZFTA  /  (F27.5.4E22.5)  ) 

5  C 0 N T  INUF _ 

PUNCH  1003 

1003  FORMAT  I  BOM***************** ********************************* ****** 

l  ****  ******  ************** ) 

6  DO  7  1=1, NPTS 

Tl(I) =-YRAR I  I) 

7  CONTINUE _ 

PUNCH  1004,  (  XBARI  I ) , YBAR I  I ) , ZBAR I  I ), XBARI  I  ),Tl( I ) »  ZRAR  I  I  )  , 

1  1*1, NPTS) 

1004  FORMAT! 9F8«  3) 

RETURN 

END 
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